Meron Pairs and Fermion Zero Modes 
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Merons, conjectured as a semiclassical mechanism for color confinement in QCD, have been de- 
scribed analytically by either infinite action configurations or an Ansatz with discontinuous action. 
We construct a smooth, finite action, stationary lattice solution corresponding to a meron pair. We 
also derive an analytical solution for the zero mode of the meron pair Ansatz, show that it has 
the qualitative behavior of the exact zero mode of the lattice solution, and propose the use of zero 
modes to identify meron gauge field configurations in stochastic evaluations of the lattice QCD path 
integral. 

PACS number(s): 12.38.Aw, 11.15.Ha, 11.15.Kc, 12.38.Gc 



Chiral symmetry breaking and confinement of color 
charge, the two essential features of low-energy QCD, 
cannot be understood perturbatively in the coupling con- 
stant. Two complementary nonperturbative approaches 
are analytical approximation using semiclassical analysis 
of the QCD path integral (analogous to the WKB ap- 
proximation in quantum mechanics) and numerical solu- 
tion using lattice gauge theory. Combining physical un- 
derstanding from the semiclassical approximation with 
systematic rigor of the lattice has provided valuable in- 
sight into hadron structure and our goal is to use this 
combined approach to study confinement. 

Most semiclassical treatments of QCD have focused 
on instantons, leading to a qualitative, and in some 
cases even quantitative, understanding of chiral symme- 
try breaking and the structure of the QCD vacuum Q . 
However, simple configurations of instantons, such as a 
dilute gas ^ or a random superposition do not con- 
fine color charges. Of alternative semiclassical mecha- 
nisms proposed for confinement, including merons 
monopoles ||, and center vortices 0], we focus on 
merons, which are solutions to the classical Yang-Mills 
field equations with one-half unit of localized topologi- 
cal charge. They are the (3+l)-dimensional nonabelian 
extension of the 't Hooft-Polyakov monopole configura- 
tion, which has been shown to lead to confinement in 
2+1 dimensions J7j. Ultimately, we seek to understand if 
merons generate area and perimeter laws characterizing 
the behavior of confined and screened color sources. 

A purely analytic study of merons has proven in- 
tractable for several reasons. Unlike instantons, no exact 
solution exists for more than two merons Q , gauge fields 
for isolated merons fall off too slowly (A ~ l/x) to super- 
pose them, and the field strength is singular. A patched 
Ansatz configuration that removes the singularities ||] 
does not satisfy the classical Yang-Mills equations, pre- 
venting even calculation of gaussian fluctuations around 
a meron pair || . 

This letter lays the groundwork for using lattice tech- 
niques to investigate the presence and role of merons in 



QCD. We construct a smooth, finite action configuration 
corresponding to a meron pair, which with its associated 
quantum fluctuations, should be present in Monte Carlo 
lattice QCD calculations. Similar work with instantons 
has shown fermionic zero modes are a powerful tool for 
identifying topological objects on the lattice. Therefore, 
we also analytically derive the zero mode of the contin- 
uum meron pair Ansatz and compare with exact numer- 
ical results. For simplicity, we restrict our discussion to 
SU(2) color without loss of generality. 

Two known solutions to the classical Yang-Mills equa- 
tions in four Euclidean dimensions are instantons and 
merons. Both have topological charge, can be interpreted 
as tunneling solutions, and can be written in the general 
form (for the covariant derivative — — iA^a a /2 
and 't Hooft symbol rj a ^) 



(1) 



with H(x 2 ) — l/(x 2 + p 2 ) for an instanton and II(x 2 ) = 
for a meron. 

Conformal symmetry of the classical Yang-Mills ac- 
tion, in particular under inversion x M — > x^/x 2 , shows 
that in addition to a meron at the origin, there is a sec- 
ond meron at infinity. These two merons can be mapped 
to arbitrary positions, which we define to be the origin 
and d^. After a gauge transformation, the gauge field for 



the two merons takes the simple form |10 

(x 
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Similar to instantons Q , a meron pair can be expressed in 
singular gauge by performing a large gauge transforma- 
tion about the mid-point of the pair, resulting in a gauge 
field that falls off faster at large distances (A ~ x~ 3 ). 
A careful treatment of the singularities shows that the 
topological charge density is JTf| 



Q(x) 



Tr 



1 



1 



F, lJ F^)=-5\x) + -5\x-d) , (3) 



16tt 2 V J 2 " v ~ y ' 2 

yielding total topological charge Q = 1, just like the in- 
stanton. 
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FIG. 1. Meron pair separated by d = V Rr regulated with 
instanton caps, each containing | topological charge. 

The gauge field Eq. (||) has infinite action density at 
the singularities x^ = {0, cL}. Hence, a finite action 
Ansatz has been suggested [[| 



x z < r 
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< Vx 2 < i? , (4) 



R < 



x 2 +R 2 ' 

Here, the singular meron fields for Va? < r and v^c 2 " > R 
are replaced by instanton caps, each containing topologi- 
cal charge \ to agree with Eq. (||) . The radii r and R are 
arbitrary and we will refer to the three separate regions 
defined in Eq. (|J) as regions I, II, and III, respectively. 
The action for this configuration is (with So = 8tt 2 / g 2 ) 



S = S, ( 1 + ^ln- 

8 r 



(5) 



which shows the divergence in the r — ► or R — > oo limit. 
Although this patching of instanton caps is continuous, 
the derivatives are not, and so the equations of motion 
are violated at the boundaries of the regions in Eq. 

Applying the same transformations used to attain 
Eq. (||) to the instanton cap solution, regions I and III 
become four-dimensional spheres each containing half an 
instanton ||. The geometry of the instanton caps is 
shown in Fig. |l| for the symmetric choice of displacement 
d = yRr along the ^-direction. Note that the action 
for the complicated geometry of Fig. [j] is still given by 
Eq. (||). We will concentrate on this case below, and 
generalization to a different d is straightforward. 

The original positions of the two merons x M = {0, d^} 
are not the centers of the spheres, nor are they the po- 
sitions of maximum action density, which occurs within 



the spheres with S n 



(48/g 2 )(R + r) 4 /d 8 at 



(Mi 



d 2 ^' 



(Mm), 



R 2 



R 2 



T 2 ^ 



(6) 



However, in the limit r — > and R — > oo holding d fixed, 
the spheres will shrink around the original points reduc- 
ing to the bare meron pair in Eq. (^). In the opposite 
limit R — > r, the radii of the spheres increase to infinity, 
leaving an instanton of size p = d. 



An instanton can therefore be interpreted as consisting 
of a meron pair. This complements the fact that instan- 
tons and antiinstantons have dipole interactions between 
each other. If there exists a regime in QCD where meron 
entropy contributes more to the free energy than the 
logarithmic potential between pairs, like the Kosterlitz- 
Thouless phase transition, instantons could break apart 
into meron pairs j^j. The intrinsic size of the instanton 
caps would then be determined by the instanton scale p. 

The Atiyah-Singer index theorem states that a fermion 
in the presence of a gauge field with topological charge 
Q, described by the equation 



pip = \ip , with ip 



IpL 



(7) 



has riR right-handed and left-handed zero modes (de- 
fined by A = 0) such that Q = nj, — tir. This can some- 
times be strengthened by a vanishing theorem |T]]] . Ap- 
plying p twice to tj} decouples the right- and left-hand 
components, and focusing on the equation for tpR, gives 



(8) 



where a a acts on the spin indices of ipn and F^ v acts on 
the color indices. For a self-dual gauge field (like an in- 
stanton), the second term in Eq. (^) is zero; and since D 2 
is a negative definite operator, there are no normalizable 
right-handed zero modes, implying Q — n^. Although 
the meron pair is not self-dual, the second term can be 
shown to be negative definite as well, leading to the same 
conclusion. 

The finite action Q — 1 gauge field Eq. ([!]) has a 
fermion zero mode given by 



TP 



(9) 



with normalization J\f and s = it72 coupling the spin in- 
dex a to the color index a (both of which can be either 
1 or 2) in a singlet configuration. The gauge field for the 
meron pair with instanton caps corresponds to 
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for regions i = I, III, 



region II, 



(10) 
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lution to the zero mode is then Eq. (||) with Af 

1/2 



The normalized so- 
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. Note that the unregulated meron 
oo) has a normalizable zero mode 



ds° 
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FIG. 2. The action density S(x) for the regulated meron 
pair, normalized by the instanton action So = 8tv 2 /g 2 , sliced 
through the center of the configuration in the direction of sep- 
aration. Shown are the analytic (solid), initial lattice (open 
circles) and cooled lattice (filled circles) configurations. The 
same is shown for the topological charge density Q(x) and 
fermion zero mode density i/j^ip(x). 

The gauge-invariant zero mode density ip'tp(x) has a 
bridge between two merons that falls off like x~ 3 in con- 
trast to the x~ 6 fall-off in all other directions. This be- 
havior can be used to identify merons when analyzing 
their zero modes on the lattice, similar to what was done 
for instantons in Ref. [fl3|| . 

As mentioned above, the patching of instanton caps to 
obtain an explicit analytic solution has unavoidable and 
unphysical discontinuities in the action density. There- 
fore, in order to study this solution further, we put the 
gauge field on a space-time lattice of spacing a in a box of 
size Lo x L\ x L2 x L3. The gauge-field degrees of freedom 
are replaced by the usual parallel transport link variable, 



U^{x) = Pexp 



-i J x+ae f A v {z)dz v . The exponenti- 
ated integral can be performed analytically within the 
instanton caps, producing arctangents. Outside of the 
caps, the integral is evaluated numerically by dividing 
each link into as many sub-links as necessary to reduce 
the 0(a 3 ) path-ordering errors below machine precision. 

We then calculate the action density S(x) for a meron 
pair with instanton caps using the improved action of 
Rcf. |0 and the topological charge density Q{x) using 
products of clovers. These (open circles) are compared 
with the patched Ansatz results (solid line) in Fig. || for 
r = 9 and d = 20.8 (in units of the lattice spacing) on a 
32 3 x 40 lattice. We use the Arnoldi method to solve for 
the zero mode of this gauge configuration and hence the 
density ip^ip(x) : which is also compared with the analytic 
result in the same figure, showing excellent agreement. 
Two important features of the patched Ansatz also ev- 
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FIG. 3. The topological charge density Q(x) on the left 
and fermion zero mode density ip^ip(x) on the right for a 
cooled meron pair with r — 9, separated by distance d = 20.8, 
17.7, and 14.1 respectively, in the (z, t)-plane. 

ident with the lattice representation are the discontinu- 
ities in the action density at the boundary of the instan- 
ton caps and vanishing of the topological charge density 
outside of the caps as given by Eq. ([5]) . The discontinuity 
in the action density is unphysical and can be smoothed 
out by using a relaxation algorithm to iteratively min- 
imize the lattice action |lf|. On a sweep through the 
lattice, referred to as a cooling step, each link is cho- 
sen to locally minimize the action density. Since a sin- 
gle instanton is already a minimum of the action, this 
algorithm would leave the instanton unchanged (for a 
suitably improved lattice action) . The result for the reg- 
ulated meron pair after ten cooling steps is represented 
in Fig. H by the filled circles, showing the discontinuities 
in the action density have already been smoothed out. 
In Fig. ||, we also plot side-by-side the cooled Q(x) and 
ijfiipix) in the (z,£)-plane for three different meron pair 
separations. As the separation vanishes, the zero mode 
goes over into the well-known instanton zero mode result. 

Like instanton-antiinstanton pairs, however, a meron 
pair is not a strict minimum of the action, since it has a 
weak attractive interaction, Eq. (^), and under repeated 
relaxation will coalesce to form an instanton. Neverthe- 
less, just as it is important to sum all the quasi-stationary 
instanton-antiinstanton configurations to obtain essen- 
tial nonperturbative physics Kg ), meron pairs may be 
expected to play an analogous role. A precise framework 
for including quasi-stationary meron pairs is to express 
the partition function as an integral over a suitable col- 
lective variable q with a constraint Q[A] as follows 



3 



Z = JdA exp{-S[A]} = Jdq Z q , (12) 

Z q = J DA exp{-S[A}-\(Q[A}-q) 2 } . (13) 

The meron pair is then a true minimum of the effective 
action in Eq. ( |l3| ) , allowing for a semiclassical treatment 
of Z q . We choose Q[A] to be the quadrupole moment 



2.5 



3z 2 



y 2 — t of the topological charge. 



The criterion for a good collective variable q of the 
system is that the gradient in the direction of q is small 
compared to the curvature associated with all the quan- 
tum fluctuations. In this case, we have an adiabatic 
limit in which relaxation of the unconstrained meron pair 
slowly evolves through a sequence of quasi-stationary so- 
lutions, each of which is close to a corresponding sta- 
tionary constrained solution. Detailed comparison of 
the quasi-stationary and constrained solutions shows this 
adiabatic limit is well satisfied. Therefore, we will present 
here the action of a meron pair as it freely coalesces into 
an instanton. To compare with the patched Ansatz, we 
need the meron separation d and radii r and R of the 
lattice configuration. We find these by first measuring 
the separation of the two maxima in the action density 
A = | Mi — Mni| (using cubic splines) and their values 
(which are the same in the symmetric case, denoted by 
Smax = 48/g 2 w 4 ), and then using Eq. (g) to give 



Aw 2 



{r,R} 



2wd 
d±A 



(14) 



In Fig. ||, we plot the total action of a meron pair as a 
function of d/r (solid line), which for the analytic case is 
given by Eq. (^) with R — d 2 jr. Also shown are cooling 
trajectories for four lattice configurations with different 
initial meron pair separations. Each case starts with a 
patched Ansatz of a given separation and the lattice ac- 
tion agrees with that of the analytic Ansatz. 

The first few cooling steps primarily decrease the ac- 
tion without changing the collective variable (which is 
now effectively d/r), as observed above in Fig. ||. Fur- 
ther cooling gradually decreases the collective variable, 
tracing out a new logarithmic curve for the total action 
(dashed line in Fig. |]), which is about 0.25 lower than the 
analytic case. This curve represents the total action of 
the smooth adiabatic or constrained lattice meron pair 
solutions. The essential property of merons that could 
allow them to dominate the path integral and confine 
color charge is the logarithmic interaction Eq. (||) which 
is weak enough to be dominated by the meron entropy {1 . 
Hence, the key physical result from Fig. [| is the fact that 
smooth adiabatic or constrained lattice meron pair solu- 
tions clearly exhibit this logarithmic behavior. 

In summary, we have found stationary meron pair so- 
lutions, shown that they have a logarithmic interaction, 
and shown that they have characteristic fermion zero 
modes localized about the individual merons. Hence the 
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FIG. 4. Action of a regulated meron pair as a function of 
d/r for initial lattice configurations of dinit = 12, 14, 18, 22 
and nnit = 4 on a 32 3 x 40 lattice. The circles show cooling 
trajectories with each cooling step marked. 

presence and role of merons in numerical evaluations of 
the QCD path integral should be investigated, and the 
study of fermionic zero modes and cooling of gauge con- 
figurations provide possible tools to do so. 
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